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ABSTRACT 



The Interstellar Boundary Explorer (IBEX) spacecraft is providing the first 
all-sky maps of the energetic neutral atoms (ENAs) produced by charge-exchange 
between interstellar neutral H° atoms and heliospheric solar wind and pickup ions 
in the heliosphere boundary regions. The 'edge' of the interstellar cloud presently 
surrounding the heliosphere extends less than 0.1 pc in the upwind direction, ter- 
minating at an unknown distance, indicating that the outer boundary conditions 
of the heliosphere could change during the lifetime of the IBEX satellite. Using 
reasonable values for future outer heliosphere boundary conditions, ENA fluxes 
are predicted for one possible source of ENAs coming from outside of the he- 
liopause. The ENA production simulations use three-dimensional MHD plasma 
models of the heliosphere that include a kinetic description of neutrals and a 
Lorentzian distribution for ions. Based on this ENA production model, it is 
then shown that the sensitivities of the IBEX 1.1 keV skymaps are sufficient to 
detect the variations in ENA fluxes that are expected to accompany the solar 
transition into the next upwind cloud. Approximately 20% of the IBEX 1.1 keV 
pixels appear capable of detecting the predicted model differences at the 3a level, 
with these pixels concentrated in the Ribbon region. Regardless of the detailed 
ENA production model, the success of the modeled B • R ~ directions in re- 
producing the Ribbon locus, together with our results, indicate that the Ribbon 
phenomenon traces the variations in the heliosphere distortion caused by the 
relative pressures of the interstellar magnetic and gaseous components. 

Subject headings: ISM: magnetic fields, clouds, HI — solar system: general — stars: 
winds, outflows 
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Introduction 



The dynamic heliosphere varies with the properties of the surrounding interstellar cloud 
and the solar wind. A theoretical study of the expected heliosphere response to different 
types of interstellar clouds show that both the overall dimensions and hydrogen filtration 
should vary substant ially with variations in th e physical properties of circumheliospheric 



interstellar material ( iMiiller et al. 



2006 



20081 ). The presence of an interstellar magnetic 



field causes heliosphere asymmetries that can diagnose the properties of the surrounding 
interstellar material, but which are partially offset by the charge-exchange coupling of 



interstellar H° and protons upstream of the h eliopause (IPogorelov et al. 



2009 



Ratkiewicz et al. 



2008 



Izmodenov 



2009 



2009c : 



Opher et al. 



interstellar gas inside of the heliosphere (e.g. 



The veloc ity discontinuity observed between 



Witte 



2004 ) and interstell ar material (ISM ) 



towa r ds the nearest stars in the upwinc 



1984 



Lallement et al. 



1995 



Wood et al. 



direction (36 Oph and a C en 



2000 



Linsky fc Wood 



Landsman et al. 



19961 ) is usually interpreted 



to indicate that the Sun is immersed in one interstellar cloud today, but will enter a separate 
cloud sometime in the next ~ 4000 years. The Interstellar Boundary Explorer (IBEX) is 
for the first time mapping heliospheric energetic neutral atoms, formed b y charge-exchang e 



between solar wind ions and i 


Dickup ions wi 


h neu 


;ral interstellar a 


;oms 


2009b 


a; 


Funsten et al. 


2009b 


Fuselier et al. 


2009; 


Schwadron et al. 


2009) 



(IMcComas et al. 



20091 ). In this paper we 



show that, depending on the source of EN As observed, IBEX is capable of detecting the 
variable heliosphere boundary conditions that might accompany the expected (someday, it 
could be soon) solar transition into a new interstellar environment in the upwind direction. 

The discovery of an unexpected 'Ribbon' of ENA emission, in directions where the 
interstellar magnetic field draping over the heliosphere is thought to be perpendicular to 
the sightline, showed that IBEX may be detecting plasma-neutral interactions beyond the 
heliopause. The similar spectra of EN As in the Ribbon and adjacent sightlines suggest 
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that the Ribbon represe nts a selection effect rat her t han an ENA popu 



fully independent origin. 



McComas et al 



f l2009ah and 



Schwadron et al 



at ion with an 



(12009h noted that 



the ribbon is organized by the most likely direction of the external interstellar magnetic 
field (ISMF), and proposed several different potential sources of the Ribbon including the 
possibility that the Ribbon might be created from a population of anisotropic suprathermal 
ions gyrating around the interstellar magnetic field just outside the heliopause. These 
ions could be indigenous to the outer heliosheath (beyond the heliopause) or more likely 
would arise from EN As that propagated out from the supersonic solar wind and/or inner 
heliosheath (between the termination shock and the heliopause); these authors noted that 
the problem with this idea is that the relatively confined pitch angle distributions would 
need to be maintained long enou gh to create "secondary E NAs" , which likely takes several 



years on average. Subsequently, 



Heerikhuisen et al. 



finii i 



incorporated this idea, where 



outward propagating ENAs create an anisotropic population of pickup ions (PUI) in the 
outer heliosheath that seed "secondary ENA" production several hundred AU upwind of 
the heliopause, into a quantitative model. While it is still uncertain how long the ion 
ring beam takes to scatter into a shell distribution (Florinski, Zank, Heerikhuisen, Hu and 
Khazanov, submitted), with marginally stable ring distributions predicted by some models 
for the distribution of the pitch angles of pickup ions in the outer heliosheath (Gamayhunov, 
Zhang, Rassoul, submitted to ApJ), the Heerikhuisen et al. simulation assumes that the 
re-neutralization time is essentially instantaneous compared to the scattering time so 
that the new secondary ENA will have a preferred direction that is perpendicular to the 
local ISMF direction. In this model, IBEX then sees these secondary ENAs where the 
gyration plane of the ion is aligned with the sightline to IBEX, i.e. where the sightline is 
perpendicular to the ISMF direction. 

IBEX data are uniquely qualified to simultaneously test both the direction of the ISMF 
at the heliosphere and the density of interstellar neutrals. The ISMF drapes around the 
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heliosphere, rotating by ~ 30° between 'infinity' and the heliopause near the nose. The 
Ribbon is ~ 20° wide and at least 270° long, possibly forming a complete circle in the sky. 
IBEX only sees ENAs with momentum vectors directed back towards the inner heliosphere. 
The long mean free paths of ENAs, e.g. ~ 200 AU for 1.1 keV ENAs in n(p + )~ 0.1 cm -3 
plasma, allow detection in the inner heliosphere of secondary ENAs formed in regions 
beyond the heliopause with elevated interstellar densities and relatively isotropic ENA 
velocities (compared to the outwards radial flows for primary ENAs produced in the 
supersonic solar wind, although not compared to the inner heliosheath ion populations). 
The solar wind contributing to ENA production includes both core ions from the expanding 
solar corona, and pickup ions formed inside of the heliosphere by charge exchange between 
interstellar neutrals and the core solar wind. 



The predictive capabilities of global heliosphere models have improved significantly to 
accommodate observational constraints placed by the 10 AU diffe rence in the termination 



shock distances found by the Voyager 1 and Voyager 2 spacecraft (jStone 



2007 



Opher et al 



20091 ). the ~ 5° offset between the upwind directions of interstellar H° an d He° flowing in to 



the heliosphere determined fro m SOHO/SWAN and Ulysses/GAS data (ILallement et al. 



2005 



Witte 



2004 



Frisch 



2008 



where the He° upwind direction must first be converted 



to J2000 co ordinates for this com parison), the properties of the ISM surrounding the 



heliosphere (ISlavin fc Frisch 



20081 ). and now the IBEX data on ENA fluxes and the Ribbon. 



Although the IBEX Ribbon itself was not predicted by models of ENA formation in the 
heliosphere, the global heliosphere models provide d the ISMF orientation that matches we ll 



with t he configuration of the Ribbon in the sky (ISchwadron et al. 



2009 



2009d). 



Pogorelov et al. 



In the discuss i ons be l ow we rely on t 



Heerikhuisen et al 



bQlOj); 



r e ENA production models quantified by 



Pogorelov et al. 



fl2009bf ) to predict the ENA fluxes for a 
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heliosphere immersed in the next interstellar cloud versus the present-day surrounding 
interstellar cloud. The ENA production simulations use three-dimensional MHD plasma 
models of the heliosphere that include a kinetic description of neutrals and a Lorentzian 
distribution for solar wind protons to approximate the suprathermal population of pick-up 
ions in the heliosheath region. The interstellar neutrals, that are coupled self-consistently to 
the plasma component in the MHD heliosphere models, act to symmetrize the heliosphere. 
Any asymmetry in the quiescent plasma distribution (e.g. created by ISMF) results in 
variations in the number of charge-exchange events, creating ne w ions that are decelera ted 
by the heliopause and therefore that mitigate the asymmetry (IPogorelov et al. 



2009c 



2008a. 



bj). The LISM flow, for the assumed ISMF strength, is subfast magnetosonic (Table 
1). This results in the absence of a bow shock in front of the heliopause, and increases 
the width of the region where the ISMF deviates from its unperturbed orientation. 
Since the local ISMF direction varies as interstellar protons approach the heliopause, the 
mean-free-path for the charge exchange interaction must also be included self-consistently 
in any Ribbon production models. 



The Ribbon ENAs are formed upwind of the heliopause in the 



Heerikhuisen et al. 



(120101 ) model, so the outer boundary conditions set by the ISM have a direct impact on 



predicted ENA fluxes and provide a means of estimating ENA variations from a cloud 
transition, regardless of whether the model is correct in detail. Any model that reproduces 
the observed B ■ R ~ of the Ribbon, which is seen where the ISMF (B) is perpendicular 
to the sightline (R), should provide useful insights into the deformation of the heliosphere 
due to altered boundary conditions from the v ariable interstellar wind. Our conclusions 
here rely explicitly on the assumption that the 



Heerikhuisen et al. 



( 120101 ) model provides 



a viable description of ENA production, both for the cloud we are in today and for the 
nearby cloud observed towards a Cen and 36 Oph in the upwind direction. 
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We focus on the 1.1 keV data, because the contrast between Ribbon ENA fluxes and 
diffuse ENA fluxes is stronger at this energy, and partly because of the enhanced outwards 
flow of ENAs from core solar wind ions that have typical energies near 1 keV. The ENA 
spectra are explicitly predicted by the heliosphere models, however this spectral information 
is not used here. The observed fluxes of ~ 4.5 keV ENAs are an order of magnitude lower 
than the 1.1 keV fluxes and the mean-free-paths are ~ 50% larger. Since the outflowing ion 
fluxes will decrease as ~ R~ 2 with distance from the Sun R, both parent ion densities and 
the resulting 1-AU ENA densities are lower when production regions are further from the 
Sun. 



2. Properties of the Upwind ISM 



The boundary conditions of the heliosphere are set by the ISM, and vary over 
geologically short time-scales. Interstellar clouds within ~ 50 pc flow past the Sun 
with heliocentric velocities that cluster around ~ 28 km s" 1 (after projection effects are 



removed), or ~ 1 parsec per 35,000 yrs (jFrisch fc Slavin 



20061 ). If nearby I SM is in pressur e 



equi 



ibrium, then models of the interstellar cloud around the heliosphere ( ISlavin fc Frisch 
20081 ) combined with ISM data (IRedfield fc Linskyi 120041 ) yield an estimate for the typical 
cloud length of ~ 1 pc. Cloud column densities for 23 cloud components towards stars 
within 10 pc yield a range of cloud lengths 0.06-3 pc, giving typical cloud crossing times for 
the Sun of ~ 1, 450 - 2.8 x 10 5 yearst. 

The velocity of the ISM at the heliosphere, which we term here the heliospheric ISM but 



1 These estimates assume an equilibrium thermal pressure of ~ 3 x 10~ 13 dynes cm -2 for 
the present-day cloud, uniform magnetic and cosmic ray pressures, D/H~ 1.6 x 10~ 5 , and a 
uniform proton density of 0.08 cm -3 
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which is also known as either the Local Interstellar Cloud (LIC) or the circumheliospheric 
ISM, is best set by t he velocity o f interstellar He° observed inside of the heliosphere by 



the Ulysses satellite (jWitte 



heliosheath regions (< 2%. 



200 4j) . Interstellar He° experiences minimal filtration in the 



Miiller fc Zank 



20041 ) and is not subject to radiation pressure, so 



that variation in the interstellar He° velocity as it traverses the heliosphere is minimal. The 
first comparisons between the velocities of interstellar gas in the heliosphere and outside of 
the heliosphere, toward s stars in the upwind direction, showed differences of over 3 km s _1 



( Adams fc Frisc 



19771 ). When the 26.3 km s 1 He° velocity is projected towards the 



nearest star in the upwind direction a CMa, 1.3 pc away and 50° from the heliosphere nose 
the projected heliospheric ISM velocity is -17.0 km s _1 , in contrast to th e observed cloud 
veloc ity from the unsaturated Fe + and Mg + lines of —18.0 ± 0.2 km s -1 (ILinsky fc Woo 



19961 ). The heliosphere nose direction is given by the upwind direction of He° flowing 



Wittc 



2004 . 



through the heliosphere, or \,0 = 255.4° ± 0.5°, +5.1 ± 0.2° (epoch J2000 
and private communication). The star 36 Oph is 6 pc away and 10° from the heliospheric 
nose. The projected heliospheric ISM velocity in this direction is -25.9 km s " 1 , versus th e 



obser ved cloud velocity from the Fe + and Mg + lines of —28.1 ± 0.2 km s 1 (IWood et al. 



20001 ) . The limit on a cloud component at the heliospheric ISM velocity towards 36 Oph is 



~ 6 x 10 16 cm 2 , giving an upper limit to the heliospheric ISM edge^l in this direction of 0.1 
pc which will be traversed in less than 4000 years. 

The interstellar cloud observed towards 36 Oph and a Oph is kn own as the G-cloud 



(ILallement fc Bertin 



1992 



Frisch et al. 



2002 



Redfield fc Linsky 



20081 ). since it is close to 



the Sun in the galactic hemisphere. Although other possible clou ds have been suggested as 



the next cloud to be encountered by the heliosphere ((Frisch 



20031 ). the G-cloud remains the 



2 We use n(H°)=0.19 cm 2 for the heliospheric ISM H° density, based on Model 26 in the 
radiative transfer models of ISlavin fc Frisch! ( 120081 ) . 
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most likely future heliospheric environment. 

The scenario examined here assumes that the Sun transitions directly from the 
heliospheric ISM to the G-cloud seen towards 36 Oph and a Cen. We determine the 
properties of the next cloud by assuming that the thermal and magnetic pressures in 
the heliospheric ISM are equal, and that the heliospheric ISM and G-clouds are in 
pressure equilibrium. The G-cloud temperature is found from the mass- dependent Doppler 
broadened widths of interstellar absorption lines, and is 5400 ± 500 K towards a Cen, and 
5900 ± 500 towards 36 Oph (where the cloud column density is also larger by 70%). The 
cooler G-cloud temperatures are thus compensated by neutral densities that are 20% larger 
than the heliospheric ISM. The 9% higher velocity of the G-cloud, in comparison to the 
velocity of the ISM now surrounding the Sun, will increase the interstellar ram pressure even 
if thermal pressure and ionization levels remain constant. We test the sensitivity of ENA 
emission to the ISMF direction using two separate ISMF directions for the G-cloud. The 
first assumption is that the directions of the ISMF in the heliospheric ISM and G-clouds 
are the same (Models 2 and 3 in Table 1). For Model 4 we make an arbitrary^ assumption 
for the G-cloud ISMF direction, which is 28° different from today's field but less than the 
~ 30° rotation of the ISMF between the ISM and the heliopause for the upwind direction. 
The detailed heliosphere boundary conditions used here for the next-cloud are listed in 
Table 1. We show below that the IBEX Ribbon is highly sensitive to even small variations 
in the ISMF direction, even when increased H° densities mitigate the influence of the ISMF 
on heliospheric asymmetries. 



3 Caveat: This direction, which is arbitrary from an interstellar viewpoint, was selected to 
lie on the great circle that divides the hot and cold poles of the cosmic m icrowave bac kground 
dipole moment, and which passes through the heliosphere nose region (IFrischl 120101 ). 
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3. ENA fluxes from encounter with Next Interstellar Cloud 

The heliosphere model has been run for the interstellar boundary conditions listed 
for Models 1-4 in Table 1. The 1.1 keV ENA fluxes predicted by the resulting models 
are displayed in Fig. CD Model 1 corre sponds to the heliospher e model displayed in 



Schwadron et al. 



(120091 ) . Model 2 (from iHeerikhuisen et al.ll2010l ) is the updated model we 



use for the heliosphere today. Models 3 and 4 correspond to the anticipated heliosphere 
environment in the next interstellar cloud, with Model 3 having a similar ISMF configuration 
as Model 2, and Model 4 showing a ISMF field with a different direction. The increased H° 
density in Model 3 produces a brighter ribbon, while the ram and thermal pressures slightly 
increase and the magnetic pressure slightly decreases. The magnetic field in Model 2 plays 
a bigger role in deforming the heliosphere than in Model 3. Model 4 has the Ribbon shifted 
significantly because of the different ISMF direction. 

The differences between the predicted ENA fluxes for Model 2 (the 'today-cloud') 
and Model 3 (the G-cloud assuming the same ISMF direction as today-cloud) are directly 
displayed in Fig. [2j The ENA flux differences are obtained by subtracting the predicted 
fluxes of Model 2 (F Mod2 ) from the predicted fluxes of Model 3 (F Mo d3, left). The most 
significant difference between the two models is the ram pressure of the neutrals, which is 
a factor of 1.8 larger for Model 3 versus Model 2. Over most of the sky, the higher flux 
of interstellar H° into the heliosphere for Model 3 generates larger ENA fluxes, with the 
differences approaching the brightest observed fluxes. However the red pixels in Fig. 2, left, 
show regions where the today-cloud has higher fluxes than the next-cloud, and represent the 
small shift in the Ribbon position due to the increase in the ratio of thermal to magnetic 
pressures, Pthcr/-Pma g , in Model 3 compared to Model 2. The effect of increased H° densities 
and thermal pressures are also seen in the increased ENA emissivity of the eastern flank of 
the heliosphere in the nose direction, where there is a bulge in total pressure (magnetic and 
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thermal, see Fig. 1 in iHeerikhuisen et al.ll2010f ). Fig. [21 right, shows the flux differences 



between Model 2 and Model 4, where the ISMF direction in the G-cloud has been allowed 
to vary by 20°. This modest variation in the ISMF direction, while retaining the same field 
strength, leads to an obvious shift in the Ribbon location, and shows that variations in the 
direction of the ISMF draping over the heliosphere should be apparent in the ENA data. 

The ENA-production model used here predicts that the ribbon and non-ribbon regions 
respond differently to variations in the interstellar density, because the ribbon also directly 
traces heliosphere asymmetries created by the ISMF. Fig. [3] shows the model differences 
modified by the ratio of n(H°) in Models 2 and 3 (0.65). The left figure shows 0.65 *-FMod3~ 
-^Mod2, and the right figure shows 0.65*FM O d4~ Fyioda- The background blue regions, where 
difference counts are ~ 0, indicate regions where the ENA fluxes are linearly related to the 
neutral interstellar density. The variations in the ribbon colors show that the ribbon ENAs 
trace the distortion of the heliosphere, which depends on the asymmetries introduced by 
the relative interstellar, ram, and thermal pressures, and which changes with different ISM 
conditions. 

The differences in the ENA fluxes predicted by the two next-cloud models are quite 
obvious in the high-flux regions of the Ribbon, but less obvious (for this color scale) for the 
directions towards the tail. In order to emphasize the differences in the weaker diffuse ENA 
emission originating in the low-flux tail region, Fig. H] shows the percentage differences 
between Models 3 and 2, (i 7 Mod3--^Mod2)/-^Mod2, and Models 4 and 2, (F Mod 4-F Mod2 )/FMod2, 
with enhanced color scales. The percentage differences in the tail region for Model 3 are 
larger than for Model 4. The outer heliosheath region around the tail is the most disturbed 
part of the model results, since the flows are subsonic. Hence larger relative variations in 
ENA fluxes are possible. 

The capability of IBEX to detect the ENA variations shown in Figs. [2H3] rests on 
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the predicted differences in the modeled ENA ffuxes compared to the flux uncertainties 
for the IBEX data. For this comparison, we use the first 1.1 keV ENA flux maps in the 



third en ergy passbanc 



imager (IFunsten et al. 



,- ESA 3") of the six energy passbands in the IBEX-HI neutral atom 



l2009af ). where fluxes are an order of magnitude larger than at 4.5 keV 
(ESA 6). Fig. [SI left, shows ESA3 fluxes, after correction for the Compton-Getting (CG) 
shift in the energy and spatial di stribution of high veloci ty particles due to the 30 km s _1 



orbital motion of the Earth (e.g. 



Gleeson fc Axford 



19681) H IBEX pixels are ~ 7°. The CG 



corrections are based on a power-law energy spectra that are derived from adjacent energy 
passbands, typically ~ E~ ls , which is evaluated over the look-direction and convolved with 
the energy response of the detector (see the Appendix in McComas et al., 2010, submitted 
to GRL, for details on the CG correction to the ENA energies measured by IBEX). The 
la uncertainties (dFi a ) on these fluxes have been determined from the Poisson statistics 
propagated from the measurement uncertainties. The IBEX data are built from data 
processed by the IBEX Science Operations Center (ISOC) for the first all-sky IBEX map 
(orbits 11-34). The uncertainties in the ESA3 fluxes are shown in Fig. [5j right. 

The measurement uncertainties for 1.1 keV fluxes can be compared to the predicted 
ENA flux differences for the next-cloud v ersus the present cloud, e.g. |AF mo d e i| = |-^Mod3 _ 



Mod2| 



based on the 



Heerikhuisen et al. 



(120101 ) model. For the comparison we preselect 



data points with signal-to-noise S/N> 3. In Fig. [61 individual pixels in the Ribbon 
region for both models have values of | AF mo d e i|/di 7 i ->> 5. The same difference map is 
plotted in Fig. [7J but now color-coded to enhance the differences in the tail region. Lower 
ENA fluxes towards the tail yield | AF moc iei|/(i-Fi -~ 1 — 2 for individual pixels, which is 
somewhat larger for Model 3 than for Model 4. Groups of 25 pixels would yield a factor of 



4 We have used the IBEX Compton-Getting corrected data set "flxset_hd60-id-base-0071- 
2010-04-09. sav", that is available at the IBEX Science Operations Center (ISOC). 
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5 improvement in the S/N of the difference maps, while effectively smoothing the data over 
~ 125 square-degrees, and still should provide a significant test. In order to use ENAs from 
the tail for identifying the next-cloud, either pixels in the tail must be grouped to improve 
statistics, or the comparison should wait for the better statistics of future skymaps. 

The predicted ENA flux differences between the today-cloud and next cloud are 
testable with IBEX data. Twenty percent of the ESA 3(1.1 keV) pixels with signal-to-noise 
S/N> 3 test the flux differences between Model 3 and Model 2 at the 3a level, or 
|AF mo dei|/<ii 71 icr> 3. In addition, 49% of the pixels test these flux differences at the la level, 
with |AF modc i|/(iF 1(7 > 1 (Fig. [8]). Similar values are found for comparisons between the 
predicted flux differences between Model 4 and Model 2, where 18% of the pixels show 
model differences that are larger than the 3a ESA3 flux uncertainties. 

We have also evaluated the variations in the 4.5 keV ENA fluxes for the environment 
of the next cloud (Fig. [9]). Although the count rates in IBEX-HI ESA 6, at 4.5 keV, are 
lower by an order of magnitude than at 1.1 keV (Fig. UM left), flux variations are predicted 
to occur when the effect of the increased interstellar density and velocity are included 
(Model 3 vrs. Model 2), as well as when the ISMF direction varies (Model 4 vrs. Model 
2). For example, Model 2 has fluxes 5-6 counts cm -2 sec -1 sr _1 keV -1 at the locations 
A,/3=253°, — 33° and A,/3=268°, 8°. For the same locations, Model 3 has fluxes a factor 
of ~ 2 higher. Further study of the energy spectrum of heliosheath ions is in progress, 
however, since the relatively large tail brightness predicted at 4.5 keV by Model 2 is difficult 
to distinguish in the data. 
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4. Discussion 

If the Ribbon ENAs are produced as secondary ENAs beyond the heliopause as 
suggested by McComas et al. (2009) and Schwadron et al. (2009), and quantified by 
Heerikhuisen et al. (2010), then the position and intensity of the Ribbon provides a robust 
diagnostic of the interstellar magnetic field direction, neutral densities, and the cloud 
ram pressure at the heliosphere. Differences between the velocity of ISM inside of the 
heliosphere and towards the nearest stars in the upwind direction indicate that the Sun 
is near or at the edge of an interstellar cloud. Based on observations of the upwind ISM, 
and assuming that the upwind cloud is in pressure equilibrium with the heliospheric ISM, 
the next cloud is modeled with densities that are ~ 50% larger, and a heliocentric velocity 
9% larger, than the cloud today. We predict the flux of ENAs from the next interstellar 
cloud to surround the heliosphere, and compare those predictions with the measurement 
uncertainties of the 1.1 keV ENA fluxes detected by IB EX in its first 6 month sk ymap. 
These results rely explicitly on the assumption that the llieerikhuisen et al. f 201o| l model 
is a viable description of ENA production both for the cloud we are in today and for the 
nearby cloud observed towards a Cen and 36 Oph in the upwind directio n. Although our 



Heerikhuisen et al. 



(120 Iff ) model, this study 



detailed conclusions rely on the accuracy of the 
is a useful gedanken experiment that will help us understand the Ribbon sensitivity to 
variations in the properties of the ISM around the heliosphere. 

The variations in ENA fluxes predicted for entry into the next cloud significantly exceed 
the measurement uncertainties for 20% of the ESA 3 pixels, which tend to be concentrated 
in the upwind hemisphere, and the variations are larger by a factor of two in some regions. 
The variations occur because the relative contributions of magnetic pressure and thermal 
ram pressure that deform the heliosphere are sensitive functions of the boundary conditions 
imposed by the ISM. If, in addition, the direction of the interstellar magnetic field shifts by 
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as much as 28°, which is slightly larger than the Ribbon width, then significant differences 
in the ENA fluxes should be observed in individual pixels near the Ribbon. The heliosphere 
regions with very low ENA fluxes, such as the tail, provide a test of the cloud properties 
only if pixels are combined for better statistical significance before comparison with model 
predictions. The long term capability of IBEX to realistically detect such variations, which 
will also be superimposed on possible solar cycle variations, requires that the efficiencies in 
the IBEX sensors (conversion, scattering, sputtering in the conversion subsystem, secondary 
electron emission at the detector foils, microchannel plate efficiencies, for instance) either 
remain stable over years, or alternatively that the instrument performances are tightly 
monitored. IBEX-HI detector sensitivity is continuously monitor ed in a number of w ays, 



20051 ). and 



such as comparison of coincident and non-coincident count rates (IFunsten et al. 
periodic gain tests. 

Variations in the energy dependence and fluxes of ENAs will occur because of the 
variation of solar wind properties over the solar cycle. These solar cycle contributions 
fortunately can be modeled in detail using past and present data on the solar wind, and 
models of the heliosphere response to these variations. Every IBEX skymap is a historical 
map of the solar cycle because of the energy dependence of ENA travel times and cross 
sections (McComas et al. submitted, 2010), so that unraveling the solar cycle dependence 
will simultaneously constrain the ENA production models and improve future predictions 
of the ENA variations due to the next interstellar cloud. 

In this discussion we considered the scenario where the next cloud is faster, slightly 
cooler, and more dense than the heliospheric ISM gas, as expected from pressure equilibrium 
and observational data. This study is a proof-of-concept, since the properties of the cloud 
edges are not established. The more extreme possibilities for the next galactic environment 
to be encountered by the Sun include a hot plasma without neutrals, and a cloud interface 
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that is either evaporative or mixed with hot plasma by shear flows. If interstellar clouds 
within 10 pc are in pressure equilibrium, they will typically fill 20% of the sightline to 
the stars. The intervening voids evidently will be filled with the low density hot gas that 
creates the Local Bubble soft X-ray emission, although t he emissivity of this plas ma is 



somewhat uncertain because of solar wind contamination ( iKoutroumpa et al 



20081). Should 



the heliosphere enter the diffuse plasma attributed to the Local Bubble interior, both 
interstellar neutrals and exo-heliospheric ENAs will vanish. IBEX and other spacecraft will 
readily detect this condition. Another possibility for the next solar environment would be 
an evaporative interface that would form upwind between the heliospheric ISM and hot 
plasma. Such an interface will show steep increases in the cloud velocity and pressure, and 
decreases in densit y, over spatial scales t hat are determined by the angle between the ISMF 



and cloud surface (jSlavin fc Frisch 



2008 . 



The present study considers the ENA fluxes for two separate models of the 
circumheliospheric cloud, but ignores possible variations due to changes in the heliosphere 
configuration during the transitions between the two clouds. The predicted thickness of the 
conductive boundary on the cloud around the heliosphere, defined as where the temperature 
falls to 50% of the asymptotic temperature, i s 0.32-0.34 pc for an ISMF direction that 



makes an angle of 30° with the cloud surface (jSlavin fc Frisch 



2008 



Models 26 and 27, also 



see Fig. 2, where the cloud edge starts at 3 pc). In the upwind direction the outflow speeds 
in the conductive boundary are 20-30 km s _1 , and opposite to the cloud motion, so that 
the Sun could traverse the conductive boundary in approximately ~ 12,000 years for these 
models. Based on the above models, we expect the change in heliosphere properties for such 
an environment to be clearly observable in the resulting ENA flux detected by IBEX. A 
turbulent mixing layer will also produce st rong gradients in the temperature and ionization 



of the surrounding ISM (jSlavin et al 



19931 1 . The ENA emissions for a cond uctive boundary 



on the surrounding cloud are discussed in detail in 



Grzedzielski et al. 



(120101 ). where the Sun 
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is estimated to emerge from the interface within ~ 500 years. An alternative possibility is 
that the G-cloud may be denser than has been assumed in this study. If the interstellar 
N(Ca + ) absorption formed at the G-cloud velocity is entirely within a few parsecs of the 
Sun, then comparisons between the clouds in the a Ce n and a Oph sightlines suggest a tiny 



cold cloud in addition to the warmer gas (jFrisch 



2003 1 



The comparisons in this paper are made without consideration of the solar cycle, 



pressure and magnetic field that characterize the solar activity cycle (IWashimi & Tanaka 


1999: 


Scherer & Fahr 


2003; 


Zank & Muller 


2003; 


Poeorelov et al. 


2009a 


). Although the 



solar cycle will cause the heliosphere to expand and contract as the solar wind dynamic 
pressure changes, these pulses travel only a relatively short distance upstream of the 
heliopause (~ 100 AU). The influence of the solar cycle on ENA production and the Ribbon 
phenomenon is not yet understood. The Ribbon intensity may vary over latitudes due to 
the ion energy differences and travel times. The extremely low levels of solar activity during 
the first year of IBEX observations suggests the solar activity cycle variations must first 



be understood before reaching a cone 



( jPogorelov et al. 



2008b 



Sternal et al. 



usion that we have entered a new interstellar cloud 



20081 ). As the theoretical models of ENA production 
become increasingly robust, we expect that studies such as this will yield definitive 
information on both the heliosphere boundary conditions and the physical properties of the 
interstellar cloud around the Sun. Finally, while this study has examined only one of the 
possible sou rces of the Ribbon cur rently under discussion, the other ideas for producing 



the ribbon flMcComas et al. 



2009al ) also generally invoke and seek to match up with the 
orientation of the external IMF, so even if another explanation eventually becomes accepted, 
it may still be possible to directly detect the interstellar transition with IBEX. 
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Table 1: Properties of Circumheliospheric Interstellar Material used for Models 



Quantity 


LIC b 


Model l c 


Model 2 d 


Model 3 e 


Model 4 f 






Original 


Today 


Next (same ISMF) 


Next (new ISMF) 


n(H°) (cm- 3 ) 


0.19 


0.15 


0.15 


0.23 


0.23 


n(p + ) (cm -3 ) 


0.06 


0.06 


0.06 


0.08 


0.08 


V (km s -1 ) 


-26.3 


-26.4 


-26.4 


-28.8 


-28.8 


T(K) 


6300 


6530 


6530 


5400 


5400 


|£| (/iG) 


[2.7] g 


3. 


3. 


2.8 


2.8 


5 direction, A,/3 




237.°, 30° 


224°, 41° 


224°, 41° 


252°, 42° 


B direction, £,b 




22.°, 41° 


36°, 53° 


36°, 53° 


40.°, 32° 


Magnetosonic Mach 


1.0 


1.1 


1.1 


0.8 


0.8 



number 



a The 1 AU solar wind parameters used for all models arc rt(p + )=7.4 cm 3 , T=51,100 K, V=450 km s 1 , 
-Bradiail = 37.5 fiG. The ISMF direction indicated by the center of the Ribbon corresponds to A, (3 = 
221°, 39° (Funsten et al. 2009). 

b These values f o r the IS M form i ng th e heliosphere boundary conditions are based on Model 26 in 



Slavin fc Frischl |2008l) and 



Schwadron et al 



Wittd (2004) 



d This model reproduces the IBEX ribbon ([Heerikhuisen et al 



(|2009l ) used this model (|Pogorelov et al.ll2008al ) in the initial analysis of IBEX data 



2010) 



e The next-cloud model, assuming the same ISMF direction as the today-model, Model 2. 
^The next-cloud model, assuming an ISMF direction that differs from Model 2. 
^Determined by assuming that thermal and magnetic pressures are equal. 
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Fig. 1. — Four models of ENA fluxes at 1.1 keV, based on different values for the ISM that 



constrains the heliosphere (see 



IBEX results ( Schwadron et al. 



able 1). Model 1 is the initial model used to evaluate the 



20091 ) . Model 2 is the assumed benchmark model for the 
production of ENAs observed by IBEX today (Heerikhuisen et al. 2010). Models 3 and 4 
utilize the same model, but with different heliosphere boundary conditions appropriate for 
the next upwind cloud. Model 3 represents 1.1 keV ENAs for a heliosphere constrained by 
the physical parameters of the next cloud in the upwind direction, where it is assumed that 
cloud is in pressure equilibrium with the circumheliospheric gas. The same ISMF direction 
of A,/3 = 224°, 41° is assumed for Models 2 and 3. Model 4 is the same as Model 3, except 
that the direction of the ISMF differs by 28°, and is directed towards \,p = 252°, 42°. The 
dashed lines, in this and subsequent figures, intersect at the longitude of the heliosphere nose 



(A = 255.4°) and the ecliptic plane. 
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Model 3 - Model 2, 1 keV Model 4 - Model 2, 1 keV 

-210 60 340 610 880 -290 -30 240 500 760 




300 200 100 300 200 100 

Ecliptic Longitude Ecliptic Longitude 



Fig. 2. — The differences between the ENA fluxes produced by the nominal G-cloud and 
the nominal LIC. Left: Model 3-Model 2 (F Mo d3- F Mod2 ). Right: Model 4-Model 2 (F Mo d4- 
^Mod2)- The different distortions of the heliosphere caused by variations in the relative 
magnetic and thermal ram pressures of the ISM are apparent on the northeast flank of 
the heliopause in the left figure. The right figure shows that the ribbon configuration is 
highly sensitive to variations in the ISMF direction. The differences in the distortion of the 
heliosphere towards the tail region are marginally visible. Fluxes are given in the units of 
counts cm -2 sec -1 sr _1 keV -1 . 
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H-Normalized Mod 3 Flux--Mod 2 Flux, 1 keV H-Normalized Mod 4 Flux-Mod 2 Flux, 1 keV 

-280 -80 120 310 510 -350 -150 60 270 480 




300 200 100 300 200 100 

Ecliptic Longitude Ecliptic Longitude 



Fig. 3. — The differences between the ENA fluxes produced by the nominal G-cloud and 
the nominal present-day cloud, after normalizing the fluxes in Models 3 and 4 by the ra- 
tio of the neutral densities (0.65). Left: (n(H ) Mod 2/n(H ) Mod3 )*F Mod 3- F Mod2 . Right: 
(n(H°) Mod2 /ri(H ) Mod 4)*F Mod 4- F Mod2 . Zero values of the difference (blue) along the rib- 
bon imply that the ribbons in the two models overlay each other. The background blue 
regions, where difference counts are ~ 0, indicate the regions where ENA fluxes are linearly 
related to the neutral interstellar density in these models. Fluxes are given in the units of 
counts cm -2 sec -1 sr -1 keV -1 . 
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Fig. 4. — The percentage differences between the 1.1 keV ENA fluxes produced by the 
nominal G-cloud and the nominal present-day cloud. Left: 100(FM O d3^-^Mod2)/-pMod2- Right: 
100(F Mo d4 _ -pMod2)/-fMod2- The right figure shows that the ribbon configuration is highly 
sensitive to any variation in the ISMF direction. 



-30 - 




Fig. 5. — Left: The first 1.1 keV IBEX ENA skymap, corrected for the Compton-Getting 
effect (see text). Right: The la statistical uncertainties on these data (see text). 
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1 Model 3- Model 2|,'1-3fl Flu* Rallotof 1 kaVSfN>3 | Model 4 ■ Model 2|/1-sig F*w Ratio tor 1 kaV&N>3 

1 I 13 19 2* 1 fl 13 IS 2' 




300 200 100 300 2O0 100 

Ecliptic Longitude Ecliptic Longitude 



Fig. 6. — The differences between ENA ffuxes predicted by models constrained by different 
boundary conditions, and compared to the la measurement uncertainties in the ESA 3 
energy passband (1.1 keV, see text). Left: The ratio of the absolute value of Model 3 minus 
Model 2 fluxes, to la uncertainties (|i 7 Mod3^-fMod2|)/^io-)- Right: Same as left figure, but 
showing the absolute value of Model 4 minus Model 2 fluxes (\FMod4.-FMod2\) / dFi a ) . Pixels 
are left blank where fluxes are insignificant (S/N< 3). 
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Fig. 7. — Same plot as Fig. El except that the color scale is coded to enhance the low flux 
areas in the tail region. 
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| Model 3-Model 2 1 / 1 sig ESA3 | Model 4-Model 2 1 / 1 sig. 



Fig. 8. — The ESA 3 fluxes (ordinate) are compared to the significance of the model dif- 
ferences (abscissa) at each pixel on the sky for Models 3 and 2 (left) and Models 4 and 2 
(right). The abscissa shows the differences in the fluxes of the two models, divided by the 
la measurement uncertainties on the ESA 3 data. Only "good" ESA 3 pixels with S/N> 3 
are plotted. The differences between Models 3 and 2 are tested at the 3a level by 20% of all 
ESA 3 pixels, while 49% of all pixels test these differences at the la level, for example. 
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4.5 keV Fluxes (counts cm" 2 sec" 1 sr" 1 keV" 1 ) 
3 7 11 15 




300 200 100 300 200 100 

Ecliptic Longitude 



Fig. 9. — The predicted ENA fluxes at 4.5 keV (ESA 6 energy passband) are shown for 
Model 2 (left) and Model 3 (right). 
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Fig. 10. — The observed ESA 6 ENA fluxes at 4.5 keV are shown (left), together with the 
la uncertainties on those fluxes (right). 



